Realizing singlet-triplet qubits in multivalley Si quantum dots 



Dimitrie Culcer, 1 Lukasz Cywiriski, 1 Qiuzi Li, 1 Xuedong Hu, 2 ' 3 and S. Das Sarma 1 

1 Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park MD20742-4111 

2 Joint Quantum Institute, Department of Physics, 
University of Maryland, College Park MD20742-4111 
3 Department of Physics, University at Buffalo, SUNY, Buffalo, NY 14260-1500 

There has been significant progress in the implementation and manipulation of singlet-triplet 
qubits in GaAs quantum dots. Given the considerably longer spin coherence times measured in Si, 
considerable interest has been generated recently in Si quantum dots. The physics of these systems 
is considerably more complex than the physics of GaAs quantum dots owing to the presence of the 
valley degree of freedom, which constitutes the focus of this work. In this paper we investigate the 
physics of Si quantum dots and focus on the feasibility of quantum coherent singlet-triplet qubit 
experiments analogous to those performed in GaAs. This additional degree of freedom greatly 
increases the complexity of the ground state and gives rise to highly nontrivial and interesting 
physics in the processes of qubit initialization, coherent manipulation and readout. We discuss the 
operational definition of a qubit in Si-based quantum dots. We find that in the presence of valley 
degeneracy a singlet-triplet qubit cannot be constructed, whereas for large valley splitting (^> feT) 
the experiment is similar to GaAs. We show that experiments on singlet-triplet qubits analogous 
to those in GaAs would provide a method for estimating the valley coupling in Si. A Zeeman 
field distinguishes between different initialized states for any valley splitting and provides a tool to 
determine the size of this splitting. 



I. INTRODUCTION 

Spin-based qubits are seen as promising candidates 
for scalable quantum computation, with donoi^ and 
quantum-dotp] spins at the focus of research. Elec- 
trical readout and control of single spins in quantum 
dots (QDs) have proven challenging, yet GaAs dou- 
ble quantum dots (DQDs), where spin blockade^ and 
charge s ensors^ e nable observation of single/two-spin 
dynamic d 4 * 11 * 12 * 1 ^ , have seen impressive experimental 
progress. In this article we establish the precise criteria 
for realizing spin qubits in Si QDs, where the multival- 
ley structure of the ground state introduces fundamental 
complications in distinguishing spin and orbital degrees 
of freedom. 

The original proposal by Loss and DiVincenzcP made 
use of a quantum dot array, in which one electron spin 
on each dot constitutes the qubit. More recently there 
has been significant progress in implementing an alter- 
native scheme, 4 in which the singlet and triplet states 
of two electrons in a DQD make up the qubitP^ 7 * 8 * 9 * 10 * 
One particular successful experiment involves initializa- 
tion, manipulation, and measurement of two-spin singlet 
and triplet states^. Here a (0,2) singlet state is initialized, 
where (n,m) indicates the occupancy of the left and right 
dots. Since the (0,2) — single dot — singlet and triplet are 
separated by an meV gap, initialization of the singlet 
is easy and reliable. Tuning the gate voltages then al- 
lows tunneling of one electron to the left dot to form a 
(1,1) singlet. When the bias is pushed deep into the (1,1) 
regime [where (1,1) is by far the electrostatic ground state 
configuration], the singlet and triplet are essentially de- 
generate due to the small tunnel coupling between the 
dots, so that a small magnetic field inhomogeneity (e.g. 



due to the Overhauser field of the nuclei) between the 
dots can rotate the two-electron states between the sin- 
glet and the triplet. After some mixing time in the (1,1) 
regime, tuning the bias returns the system to the (0,2) 
configuration, where electrical readout is possible due 
to spin blockade^. This experiment clearly illustrates 
the existence of quantum coherence in the DQD system, 
and the distinct possibility of using the two-electron sin- 
glet and unpolarized triplet as the two states of a logical 
qubit, with reliable initialization, single-qubit rotation, 
and measurement. 

Silicon is often regarded as the best semiconducting 
host material for spin qubits because of its excellent 
spin coherence properties: spin-orbit coupling is very 
small, the hyperfine interaction can be reduced by iso- 
topic purification^, and the electron-phonon interaction 
is weak as well. Furthermore, the mature Si microfab- 
rication technology will help attempts to scale up a Si- 
based quantum computer (QC) architecture. At present, 
Si/SiG<iP and Si/SiO^ quantum dots, and SiiF^l are 
being actively investigated and progress has been made 
in spin blockade in Si quantum dotsP^ The biggest ob- 
stacle to spin QC in Si is valley degeneracy: bulk Si has 
six degenerate conduction band minima. While this de- 
generacy can be reduced by strain or the presence of 
an interface, it complicates the orbital and spin state 
spectrum^ and leads to valley- interference effects for 
spin interactions^. At the Si/Si02 interface only two val- 
leys are relevant to the ground orbital state. Scattering 
at the interface further lifts the valley degeneracy by pro- 
ducing a valley-orbit coupling A. The magnitude of A is 
generally not known a priori and is sample-dependentP^. 
Currently measurement of valley splitting A is generally 
done for 2D electron gases at high magnetic fields, and 
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the zero-field valley splitting is extrapolated^. We note 
that the case of large valley splitti ng ha s been examined 
in a number of recent publications! 23 * 24 ^ 

In this article we study the physics of Si-based quan- 
tum dots and the feasibility of experiments analogous 
to Ref. [4 in a Si/Si0 2 (or Si/SiG<P) DQD, focusing 
on the effects of the valley degree of freedom on qubit 
initialization, operation, and spin blockade within the ef- 
fective mass approximation. We identify the conditions 
required for an operational singlet-triplet qubit in Si. We 
further demonstrate that a quantum coherent experiment 
analogous to Ref. |4 may provide a direct way to esti- 
mate the valley splitting A. While our discussion focuses 
on Si/Si02 and is directly relevant to experiments on 
Si/Si02 quantum dots 22 , the findings are generally ap- 
plicable to Si quantum dots. In addition, we expect our 
findings to be at least qualitatively relevant to other sys- 
tems in which the valley degree of freedom plays an im- 
portant role, such as carbon, in which significant progress 
has been made latelyP^^ 

The outline of this paper is as follows. We introduce 
the model of the DQD in Section [TTJ We proceed to study 
the initialization process in singlet-triplet qubits n Sec- 



tion III followed by manipulation of the qubit in Section 
IV| In Section |V| we demonstrate that a quantum coher- 
ent experiment on singlet-triplet qubits can be used to 
estimate the valley splitting. Section VI is devoted is- 



sues specific to silicon, such as interface roughness and 
the need for an external inhomogeneous magnetic field. 



Finally, Section VII contains a summary of our findings. 



II. DOUBLE QUANTUM DOT 

We choose nominally z as the growth direction for the 
Si/Si02 heterostructure we consider. The two dots are 
located at Rr,l — (±^o?0,0), where R and L stand 
for right and left respectively. The Hamiltonian is H = 

H + H v , with H = (Ei=i,2 T(i) + Vq) + ^ee, where 
T is the kinetic energy operator and Vq the confinement 
potential 

V Q = (1/2) m t u?o Min[(x - X ) 2 , (x + X ) 2 ] - eEx 
+ (l/2)m t uu 2 y 2 + (1/2) m z u 2 z z 2 , 

with m t and m z respectively the in- and out-of-plane 
Si effective masses. The Coulomb interaction between 
electrons at r\ and r2 is V ee = e 2 /(e\ri — 7*2 1), where 
e = (e^ + esio 2 )/2 includes the image charge in the Si02 
layer. H v is a single-particle phenomenological coupling 
between the valleys discussed below. The electric field E 
raises the energy of the left dot with respect to the right 
dot. The confinement potential and ground state for E = 
are identical in each dot, with the single-dot potentials 
Vr,l(x) = (1/2) m t ul (xtX ) 2 . At the Si/Si0 2 interface 
the lowest valleys are at ±(k z ), with (k z ) = 0.85(2ir/asi), 
and the lattice constant as% = 5. 43 A. The ground-state 
single-electron wave functions R ZyZ and L z ^ z represent the 
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Figure 1: (a) Single dot energy levels for finite A and B such 
that 2 A > Ez. The lowest energy state is the singlet Si, fol- 
lowed by the triplet T+, the degenerate singlet S2 /triplet To 
and triplet T_ , and the singlet S3 ■ Spin orientations of elec- 
trons in triplet states are indicated by arrows, (b) Schematic 
of the biased double dot. During initialization the detuning is 
large, and the (0,2) states are lowest in energy. After loading 
the detuning is lowered and the (1,1) states are at the same 
energy as the (0,2) states - the charge transition regime. An 
inhomogeneous magnetic field mixes the singlets and triplets. 



degenerate ±(k z ) valleys on the right and left dots re- 
spectively. In the right dot (T + Vr) R z , z = £0 Rz,z, with 
Rz-z = F R (r - R R )e ±ikz ^ r - R ^u ZjZ (r - R R ), and on the 
left R — > L. The envelope functions are 



F RtL (r-R RtL ) 



^2 



^3/4( a 2 6 )l/2 



where a = 



and b = 



e 2^ e 25^ 5 
(2) 

mt u v"'«. the in ~ plane 

(Fock-Darwin radius) and growth-direction confinement 
length. The lattice-periodic Bloch function u z ^ z (r) = 
° Z K e K V w ith K reciprocal lattice vectors. The 
overlap (L ZjZ \R Z:Z ) = e~ d ^ where d = Xq/cl. Overlaps 
such as (L z \Lz), (L Z \R Z ) are suppressed by an exponen- 

tial of the form e 4 ? where Q z = ™ z —2(k z ), with n z 
an integer. Such an exponential appears in all but one of 
the matrix elements of Hq involving functions from dif- 
ferent valleys. All such intervalley terms can be neglected 
except one, discussed below. The only nonzero matrix el- 
ements of H v are (L Z:Z -\H V \L Z , Z ) = (R z ,z\Hv\Rz,z) = A, 
with A > and assuming A has the same form on 
each dot. We define also sr = (R Z ^ Z \(T + Vq)\R Zi z), 
El = (L Zi z\(T + Vo)\L ZyZ ), and the dimensionless de- 
tuning as (el — £R)/(2deo)- Diagonalizing the single- 
particle Hamiltonian with the valley coupling we obtain 
the valley eigenstates R± = (l/y/2) (R z ±R Z ) and L± = 
(l/y/2) (L z ± L z ) with corresponding energies Sr } l ± A. 
We orthogonalize these following Ref. [27 , with L± = 
(L± - gR ± )/^l-2lg + g\ with g = (1 - VT^)/^, 
for R± one swaps L R. These are the states that will 
be used henceforth. 
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III. INITIALIZATION 

We begin by studying the initialization process, which 
involves loading two electrons onto the right dot. For this 
purpose it is imperative to analyze first the spectrum of 
the doubly-occupied right dot, that is the configuration 
(0,2). The lowest-energy two-particle spatial wave func- 
tions are (/)si,S3 = and 

f (3) 
v2 

where (i) denotes the ith electron. In the ba- 
sis {4>si-> 4>S2-> 4>t •> 4>sz} the matrix elements of the 
Hamiltonian are 2sr + u + diag (— 2A, 0, 0, 2A), 

where u = / rf 3 n / d 3 r 2 R* z f R*^V ee R^R™. 
The valley-exchange Coulomb integral 

J d 3 ri j d 3 r 2 Rt^R* E ^V ee R^R^ is not suppressed by 
an exponential, however we find its value to be <C 1/ieV 
and it will therefore be assumed of no consequence 
henceforth. For A = and no external magnetic field 
£?, all levels are degenerate so that it is impossible to 
load any particular two-electron state. The spectrum 
for finite A and £?, yielding a Zeeman energy Ez with 
2 A > Ez, is shown in Fig. [T] The triplet states thus 
split into T + , To, and T_, separated in energy by Ez- 

The loading process makes use of an outside reservoir 
with Fermi energy ef which is thermally broadened by 
~ k#T. The reservoir is tuned to be on resonance with 
the lowest-energy singlet state. Here we neglect the the 
differences in tunnel couplings between the various states 
and the reservoir. Consequently the probability of load- 
ing any of the states is proportional to the Fermi distribu- 
tion at its energy. If A ^> k#T, the lowest-energy singlet 
state can be loaded exclusively. Numerically A « O.lmeV 
is sufficient to fulfil this condition at dilution refrigerator 
temperatures of T=100mK, where k#T w O.OlmeV. In 
this regime the two-electron initialization process in a Si 
DQD is identical to the GaAs DQD in Ref. g]. 

IV. QUBIT MANIPULATION 

Manipulation of the singlet-triplet qubit involves 
switching to the configuration (1,1), where singlet-triplet 
mixing is achieved through the us of an inhomogeneous 
B. We study the Hilbert space of two electrons in a Si 
DQD. We do not include the high-energy (2,0) states. 
The seven space-symmetric Hund-Mulliken (HM) wave 
functions of the lowest-energy singlets are (j>si, (j)S2, (/>S3 
and the functions 

f (4) 




Detuning 

Figure 2: DQD spectrum for a=8.2nm, b=3nm, d=2.45 and 
A=0.1meV and zero magnetic field. The top and bottom 
anticrossings each consist of two singlets (solid lines) and one 
triplet (dotted line). In the middle anticrossing each of the 
three dashed lines represents a degenerate singlet /triplet level. 



These singlet states split into three uncoupled subspaces. 
The {05 + ,053} and {</>^ - , </>si} subspaces are com- 
posed of HM wave functions where the two electrons 
are in the same valley eigenstate, while the subspace 
{^ + , <p™ ~, } consists of wave functions where elec- 
trons are in different valley eigenstates. The five anti- 
symmetric counterparts of the states in Eq. Q, denoted 
by and 0T2, are evident (clearly cj>s\ and c\s3 

do not have antisymmetric counterparts.) These triplets 
in turn split into three subspaces, with {</>J + } and {c^ - } 
single-valley HM triplets, and {</>™ + , </>™~, </>™ d } mixed- 
valley triplet states. Since the overlap between states 
from different valleys is negligible, matrix elements of the 
form (^li^ol^s ) and {(j^^H^™*) are equal. As 
a result of this the mixed singlet and triplet subspaces 
always yield the same energy eigenvalues. 

Let us examine a concrete Si DQD with a=8.2nm, 
b=3nm, d=2.45 and A=0.1meW^. The energy levels 
are plotted in Fig. [2] as a function of the dimensionless 
detuning. At low detuning there are four (0,2) high- 
energy levels, indicated by the two solid lines (represent- 
ing singlets of the form R+R+, R_RJ) and one dashed 
line (representing a degenerate singlet /triplet of the form 
R+R-). The separation of these levels is 2 A. There 
are also eight lower-energy (1,1) levels: a degenerate sin- 
glet/triplet of the form L+R+ (top solid line), a degener- 
ate singlet /triplet of the form L_R_ (bottom solid line), 
and two degenerate valley-mixing singlets and triplets of 
the form L + R_ and L_R + . At high detuning the (0,2) 
states have lower energies than the (1,1) states. As in 
Ref. [4], varying the detuning drives the energy levels 
towards an avoided crossing where (0,2) and (1,1) are 
degenerate and split by the tunnel coupling t (here the 
splitting ~ 6/xeV). 

Thus far we have assumed that the valley splitting A 
exceeds the interdot tunnel coupling t. For generality, in 
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Figure 3: DQD spectrum a=8.2nm, b=3nm, d=2.45, 
A=0.01meV, E z = 0.05 meV, and t=0.2meV. Dashed lines 
indicate Zeeman-split triplets corresponding to T± in Fig. [I] 
Other Zeeman-split levels are not shown explicitly. 



Figure [3] we have shown the two-electron spectrum of a 
Si DQD when Act. Figures. [2] and [3] together demon- 
strate that, depending on the relative size of t, A, and 
Ez, the relative position of most energy levels can differ 
significantly, so that the loading and mixing dynamics 
of the two-electron states can vary dramatically. In the 
general case one must be prepared to expect an interme- 
diate situation, in which a clear separation of the energy 
levels into three branches may not occur and some of the 
energy levels may cross. 

As mentioned before, in the limit in which A ^> /c#T, 
the one-dot singlet state Si can be loaded exclusively. If 
we now tune the bias voltage as in Ref. [4] to shift to the 
(1,1) regime, the electron state will become cj)g~ (made 
from orbital states T_ and R-). This state can then 
mix with it triplet counterpart <j)^~ if an inhomogeneous 
magnetic field is present. No mixing with any other state 
is possible due to the energy separation and the fact that 
all inter-valley matrix elements of the Hamiltonian van- 
ish. When the biased voltage is tuned back to the (0,2) 
regime, the two electrons will either return to their ini- 



tial state 



or stay in the 



triplet state and get 



spin-blocked. Essentially the two-electron dynamics is 
confined to the lowest energy manifold of Fig. [2] (the low- 
est solid curve and the lowest horizontal dotted line), in 
exact analogy to what happens in a GaAs DQD in Ref P. 
That is, pulsed manipulation and measurement can be 
done reliably for singlet-triplet spin qubits in a Si DQD. 



V. MEASURING THE VALLEY SPLITTING 



an experiment analogous to the one described in RefPl 
First, we identify three loading/mixing regimes in Figs.[l] 
and [2] If Si is loaded, the system is driven to the an- 
ticrossing at the bottom of Fig. |2j as discussed in the 
previous paragraph, when the detuning is varied. If S2, 
T + , To, and T_ are loaded, the system is driven to the 
anticrossing in the middle of Fig. [2] If S3 is loaded, the 
system is driven to the anticrossing at the top of Fig. 2j 
One may in principle load any of the six states in Fig. ^ 
so that any of the three anticrossings may be involved in 
such an experiment. However, a magnetic field does not 
mix any of the three regimes of Fig. |2j thus once a state 
is loaded the experiment can be carried out as in [4] and 
the same readout process can be used. Interestingly, if 
Si, S2, S3, or Tq is loaded, the experiment is identical 
to Ref. [4]. At high B (in which singlets mix only with 
T triplets and vice versa) these states have an average 
probability of return of 1/2. If T + or T_ is loaded the 
average probability of return will be 1 at a high mag- 
netic field since they do not mix with other states. The 
average probability of return thus depends on the load- 
ing probabilities of the individual states. By studying 
the average probability of return one may estimate A by 
sweeping a uniform applied magnetic field (different from 
the inhomogeneous field mixing the singlets and triplets) 
as follows. 

Sweeping a uniform magnetic field changes the load- 
ing probability of the T + triplet, which in turn leads to 
changes in the measurable return probability. At low 
field Ez < 2A and the two electrons predominantly load 
into the ground singlet Si . Increasing the magnetic field 
will eventually bring T + below Si, with a crossing at 
Ez = 2 A. In Fig. [4] we plot the loading and return proba- 
bilities against the magnetic field for A ^> kgT in (a) and 
A « k#T in (b). In both cases the probability of loading 
Si and To will be very close to 1/2 at the crossing point 
when Ez = 2 A. Also, the return probability increases 
dramatically close to the crossing point, and the cross- 
ing point corresponds approximately to where the return 
probability reaches the mid point between its low-field 
and high- field values: For A = 0.1 meV, the crossing is 
at B = 1.76 T while the mid point is at B = 1.63 T; for 
A = 0.01 meV, the crossing is at B = 0.176 T while the 
mid point is at B = 0.17 T. The identification of the mag- 
netic field for this mid point thus gives a reliable estimate 
of the value of valley splitting 2 A = Ez- This method 
should succeed as long as 2 A >k#T. If 2 A <C /c^T, the 
return probability will not change much as we sweep the 
magnetic field, with the increased loading of T + compen- 
sated by the reduced loading into T_. Thus the overall 
change /no- change of return probability also gives a clear 
indication of whether A is larger than /c#T or not. 



The key question for a Si DQD is whether A ^> ksT. 
However, in general A is not known, so that the two- 
electron initialization comprises some uncertainty. Below 
we explore ways to determine the valley splitting using 



VI. ISSUES SPECIFIC TO SILICON 

In GaAs the inhomogeneous magnetic field is produced 
by the hyperfine interaction 4 . In Si the hyperfine inter- 
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Figure 4: (Color online): Loading probabilities of differ- 
ent levels and total return probability as a function of the 
magnetic field B for T=100mK and (a) A=0.1meV and (b) 
A=0.01meV. 

action is smaller and singlet-triplet mixing will be about 
two orders of magnitude slower than in GaAs. Using a 
nanomagnet one can design a particular field magnitude 
and direction, enabling better control of the spin qubit. 
For example an inhomogeneous magnetic field along the 
^-direction mixes only the singlet and To, whereas a field 
along the x-direction mixes only the singlet and T± . Ev- 
idently the issues discussed in this work are insensitive 
to the origin of the inhomogeneous magnetic field as long 
as this field is present in the system. 

The length scale of surface roughness, which deter- 
mines the spatial variation of A, as compared with the 
dot size and location is not precisely known. The pro- 
posed experiments will work as long as A varies over 
a length scale larger than the DQD size, or as long as 
the change in A does not lead to change in the compo- 
sitions of the valley eigenstates. We have assumed the 
same valley-orbit coupling in both QDs, thus the same 
valley splitting and eigenstates. A change in the valley 
composition of the eigenstates could lead to intervalley 
scatterings in the (0,2) to (1,1) transition, so that control 
of electron orbital states may become intractable. At the 
same time large variation of A across the DQD will hin- 
der the effectiveness of the experiment. It is imperative 



for experimental setups to ensure firstly that the inter- 
face is as smooth as possible and secondly that the DQD 
spans an area over which the interface roughness profile 
varies as little as possible. 



VII. SUMMARY 

We have studied the feasibility of initialization and co- 
herent manipulation of singlet-triplet qubits in multival- 
ley Si DQDs, demonstrating that the valley degree of 
freedom makes the physics of Si quantum dots consider- 
ably different from that of dots made out of single- valley 
systems such as GaAs. Various experimental outcomes 
are possible depending on the value of the valley split- 
ting A. For large A (i.e. A ^> fe^T) a quantum coherent 
experiment identical to Ref. [4] is feasible. For small A a 
number of different states may be initialized, leading to 
different experimental outcomes. One interesting high- 
light of our work is that, although several singlet /triplet 
states may be initialized, in general each state can mix 
with one other state, and no more. Therefore, in princi- 
ple, once a state is loaded, operations on it can proceed in 
a similar way to the scheme implemented in GaAs dots. 
For any A, sweeping a uniform magnetic field provides a 
useful method for estimating A. In fact, one very impor- 
tant consequence of our work is the proposed new method 
for estimating the valley splitting A in Si quantum dots, 
particularly when A < /c#T. Considering the difficul- 
ties inherent in proving that a certain state belongs to 
a particular valley, and thus in identifying a particular 
energy splitting with the valley splitting, it will be im- 
portant to have as many different methods as possible to 
measure/estimate the size of the valley splitting. 
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